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JS: ABSTRACT 

We present first comparisons between Light Element Primary Process (LEPP) 
j — | abundances observed in some ultra metal poor (UMP) stars and nucleosynthesis 

calculations based on long-time hydrodynamical simulations of core-collapse su- 
pernovae and their neutrino-driven wind. UMP star observations indicate Z>38 
elements include the contributions of at least two nucleosynthesis components: 
r-process nuclei that are synthesized by rapid neutron capture in a yet unknown 
site and LEPP elements (mainly Sr, Y, Zr). We show that neutrino-driven wind 
simulations can explain the observed LEPP pattern. We explore in detail the 



sensitivity of the calculated abundances to the electron fraction, which is a key 
nucleosynthesis parameter but poorly known due to uncertainties in neutrino in- 
teractions and transport. Our results show that the observed LEPP pattern can 
be reproduced in proton- and neutron-rich winds. 

Subject headings: nuclear reactions, nucleosynthesis, abundances - - Galaxy: 
abundances — supernovae: general 
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Introduction 



Elements heavier than iron are mainly produced by the slow (s-) and rapid (r-) neutron 
capture process es. In contrast to the s-process, for wh ich the astrophysical environments 



are identified (Busso et al. 



19991 ; IStraniero et all 120061 ). the r-process site(s) remains un- 



known (jArnould et al.l 120071 ). There are some indications that more th an one component 



or site contributes to the abundances of the so called r-process elements (IWasserburg et al. 



1996 



Qian fc Wasserbu rg 2001; 



Otsuki et al. 



20061 : 



Montes et al. 



Truran et al.l 12002 : iTravaglio et al.l 12004 : lAoki et al. 



2005 



20071 ). Most of the recent progress in understanding the 



origin of elements commonly associated with th e r-process is due to observations of ultra 
metal-poor (UMP) stars (see lSneden et al.ll2008l . for a recent review). The elemental abun- 
dances observed in the atmosphere of these very old stars come from a few events. These 
stars generally present a robust pattern for "heavy" elements 56 < Z < 83, in agreement 
with the expected contributi on of the r-process to the solar system, bu t show some scatter 
for "light" elements Z < 47 djohnson fc Bolt J2OO2I bneden et al.ll2008h . This suggests that 
at least two types of primary nucleosynthesis events contribute to the abundances of those 
elements. 

The first indic ation for these two nucle osynthesis components came from meteoritic data 



on 129 I and 182 Hf (IWasserburg et al.l 119961 ). The idea of the different components was later 



extended to account for all abund ances in metal-poor stars in the phenomenological model of 
Qian fc Wasserburg! ( )200ll . 120071 ) . They argued that supernovae from low- mass progenitors 
with 8M < M < 12M (H-source in their terminology) lead to all "heavy" and some 
"light" elements, and that explosions of more massive progenitors, 12M Q < M < 25M Q 
(L-source), contribute to th e remaining light A < 130 elements. A similar conclusion was 
reached independently by iMontes et al.l (120071 ) using a slightly different pattern for the 
L-source (LEPP source in their terminology). 



The abundances of star HD 122563 (jWestin et al.ll2000l ) were used to build an L-source 
pattern. This star shows mainly elements produced in the L-source with small c ontamination 



from the H-sour ce, characterized by t h e abu ndances observed in CS 22892-052 (jSneden et al. 



20031 ). While iQian &; Wasserburg! (12007J) assumed that all material above Z = 56 in 



HD 122563 was produced by the H-source, IMontes et al.l (120071 ) used a smaller H-source con- 



tribut ion responsible for all the abundance of elements heavier than Z = 62. IMontes et al. 
( 120071 ) discussed the assumption that HD 122563 is representative of the L-source by showing 
a hint of robustness of the observed pattern of UMP stars. However, more observations of 
L-source enriched (or r-process poor) UMP stars are needed before a de finitive conclusion is 
reache d. In this paper, the LEPP (o r L-source) pa t tern w as taken from IQian fc Wasserburg 



( 120081 ) based on the abundances of iHonda et al.l (120061 ). It covers elements in the range 
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38 < Z < 47 and it is assumed to be robust. 



More recently, iQian fe Wasserburg (120081 ) updated their two- component model to in- 
clude a third component producing Fe but not Z > 38 elements. The need for an additional 
component resulted from using newer higher-resolution data that included many stars with 
[Fe/H] < —3 where the amount of Fe (produced originally only in the L-source) relative to 
Sr (produced in both L- and H- sources) was too large compared to the yields based on the 
two-component model. 

The process leading t o elements with A < 1 30 (L-source) has been called in the litera- 



ture the weak r-process (Truran fe Cowan 



flWooslev fe Hoffmanlll992 



2000) , charged-particle reaction (C PR) process 



; iFreiburghaus et al.lll999clQ.iaii fe Wasserburd 



Montes et al. 



2007 ). and Lighter 
2003). In this pa- 



Element Primary Process (LEPP) (jTravaglio et al.l 12004 ; 
per, we refer to this process as LEPP because it does not make any reference to the specific 
nuclear reactions or a strophysical environment. The term LEPP was first introduced in 
Travaglio et al.l ( 120041 ) which used a galactic chemical evolution model to search for the as- 
trophysical environments producing the elements Sr, Y, and Zr. Using their s-process model 
and standard r-process contributions, they found that non-negligible a bundances of s evera l 
isotopes ( 86 Sr, 93 Nb, 96 Mo, 100 Ru, 104 Pd, 110 Cd) were still unexplained. iMontes et all (I2OO7I) 
showe d that the LEPP elemental abundances that result from using the iTravaglio et al. 
(120041 ) model are in agreement, within the observati onal error bars, with the modified abun- 
dances of HD 122563 (Fig. 5 in IMontes et all ((2003)). The contribution of the LEPP to the 
solar abundances depends on the s-process model. If further refinements to s-process mod- 
els eliminate or change the need for a LEPP contribution to the solar system abundances, 
the apparent agreement between the "solar" LEPP and the UMP observations would have 
been coincidental. In that case the LEPP may only contribute to the abundances of a few 
metal-poor stars. 



After the initial success of IWoosley et al.l (119941 ) in reproducing observed r-process 
abundances^, core-collapse supernovae and the subsequent neutrino-driven winds became 
one of the most promising candidates for the production of r-process elements because their 
extreme explosive conditions are very close to the ones needed for the r-process (see e.g., 
Hoffman et al.lll997l ; iThompson et al.ll200ll ; lOtsuki et al.l l20Q0h . Moreover, galactic chemi- 
cal evolution models favor core-collapse supernovae, since they occur early and frequently 
enough to account for the abundances observe d in old halo stars and in the solar system 
(llshimaru fe Wanajolll999l ; llshimaru et al.ll2004l ). Although the necessary conditions to pro- 



1 Recent work of iRoberts et al.l ( 20101) suggests that the high wind entropies obtained at that time could 
have been due to problems with the equation of state. 



-4- 



duce heavy elements (A > 130) are identified ( IMeyer et al.l Il992l ) (high entropies, low 
electron fractions, and short exp ansion timescales) , these are not found in the most recent 



long-time supernova simulations (|Pruet et al.ll2006l ; lArcones et al.ll2007t IWanajo et al.ll2009 
Fischer et allboici lliudepohl et al.lboioh . 



In this paper we investigate the possibility of producing the L EPP elements in neutrino- 
driven winds, as suggested by iQian fc WasserburgJ ( 120071 . 120081 ). The nucleosynthesis cal- 
culations fSect . 12 .3p are based on the spherically symmetric simulations of core-collapse 
supernovae of lArcones et al.l (|2007| ). where the evolution of the neutrino-driven wind was 
followed during several seconds (Sect. 12.21) . Since the neutrino transport in these simulations 
is only approximate ( IScheck et al.ll2006l ). the electron fraction is expected to have some 
uncertainty as explained in Sect. 12.11 A systematic study of the impact of the wind electron 
fraction on the nucleosynthesis is included in Sect. 13.11 We show that the nucleosynthesis 
based on these simulations do not yield heavy r-process nuclei (A > 130) but that LEPP ele- 
ments can be produced in the current neutrino-driven wind models. The different conditions 
to produce the observed LEPP pattern are investigated in Sect. 13.21 Finally, we conclude 
and summarize Sect. HJ 



2. Nucleosynthesis in neutrino-driven winds 

When a supernova explodes, matter surrounding the proto-neutron sta r is heated by 
neutrinos and expands very f ast reaching sometimes even supersonic velocity ( jDuncan et al. 
19861 : iThompson et al.ll200ll ). This is known as the neutrino-driven wind and can become 
neutron or proton rich. The electron fraction Y e is determined by the uncertain neutrino 
properties (energy and luminosity) in the region where neutrinos decouple from matter. 
This matter near the neutron star consist mainly of neutrons and protons due to the high 
temperatures in this region. When a mass element expands, its temperature decreases and 
neutrons and protons recombine to form a-particles. The density decreases but as the triple- 
alpha reaction co mbined with different a capture reactions a re still occurring, heavy seed 



nuclei may form ( jWoosley fc Hoffman! 11992 



Witti et al.l 1 19941 ) . The evolution once the a- 



particles start forming heavier nuclei depends on the electron fraction and it will be discussed 
in detail in Sect. 13.11 
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2.1. Wind electron fraction 



The evolution of the nascent neutron star and the properties of neutrinos emitted from 
its surface have a direct impact on the wind parameters that are relevant for nucleosynthesis: 
entropy (S), expansion timescale (r = r/t>| T ^ .5Mev) 5 an d electron fraction. S and r depend 
on neutrino energy, neu trino luminosity, and on the evolution of the proto-neutron star as 



flQian fc Wooslevlll996h : 



S oc L-^e-^R-^M QS , (1) 

r cc L- l t' 2 R ns Mn, , (2) 

where L, e are the neutrino luminosity and mean energy, respectively, and the dependence on 
the neutron star properties enters through its radius (Rob) and mass (M ns ). The evolution of 
these quantities is mainly given by the equation of state (EoS) at nuclear densities, which is 
specially uncertain at high densities, and by the amount of matter accreted onto the proto- 
neutron star, which is higher for more massive stellar progenitors. Eqs. (JTj)-(j2j) indicate that 
different progenitors with the same neutron star evolution lead to similar wind properties. 
When neutrino luminosities and energies decrease, the entropy slightly increases favoring the 
production of heavier elements, but at the same time the expansion timescale increases which 
counteracts and disfavors their production. Therefore, during the wind phase, relatively 
small changes in the nucleosynthesis can generally be expected when the absolute values of 
neutrino properties (i.e., L Ue + L Pe and e Ue +e Pe ) are changed. After the initial fast contraction 
of the proto-neutron star, its radius and mass changes only slightly for >1 s after bounce. 
The main change to S and r is then the result of the neutrino luminosity variation as a 
function of time. 

The electron fraction is, contrary to S and r, determined by the relative variations of 
electron neutrinos and antineutrinos properties. The exact value of neutrino energies and 
luminosities depend on the accuracy of the supern ova neutrino transport calculations and on 



detai ls of neutrino interactions at high densities ( iRampp fc Jankall2000l ; iMezzacappa et al. 



200ll ). At the high temperatures (T > 10 GK) where the nucleosynthesis calculations are 
started, nuclei are photo-dissociated into neutrons and protons and their abundances are 
given by Y p = Y e and Y n = 1 — Y e , respectively. The electron fraction is determined by 
charged- current reactions between neutrinos and nucleons and their inverse reactions: 

v e + n < — > p + e , (3) 

v e + p i — > n + e + . (4) 

The evolution of the electron fraction is consequently given by: 

dY e 

j 7 A Uen l n Ap e lp Ap e plp "T A ne + 1 n , yo ) 
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where Aj are the reaction rates for the forward and backward reactions in Eqs. (jSI)-©- 
Neutrino emission reactions are negligibl e in the wind bec a use th eir rates rapidly drop with 
temperature (X pe / ne + ex T 5 ). Following iQian &: Woosleyl ( 119961 ). the equilibrium initial Y e 
in the wind is given by 

-i -l 



Y P 



K 



A^, n + A 



"e,P 



1 + 



L Pe e Pe - 2 A + 1.2A 2 /e Pe 
L Ue e Ue + 2A + 1.2Aye Ue 



(6) 



Here the neutrino absorption rates are used in a simple fo rm (IQian fc Woosleylll996l ) with- 



out considering weak magnetism and recoil corrections ( [Horowitz fc Lilll999l ). L Ve is the 



neutrino luminosity, e Uf , = (e 2 )/(e) ~ A.lkT Ve is the ratio of mean squared neutrino energy 
((g 2 )) and mean neutrino energy ((e)), and kT Ue is the neutrino temperature in MeV (and 
similarly for the antineutrinos), A denotes the neutron-proton mass difference. This simple 
expression gi yes only an estima te of the Y e value found by supernova neutrino transport 
calculations (IFischer et al.l l2010h . The neutrino two-color plot shown in Fig. [1] illustrates 
the dependence of the Y e on the neutrino energies and luminosities using Eq. (Q. It also 
shows the (anti) neutrino energies obtained 10 s after bounce in win d simulations reported 



in the literature. T he most recent long-time supernova simulations (jHudepohl et al.l 12010 



Fischer et al.ll2010f ) obtain lower antineutrino energies compared to the earlier models, lead- 
ing to proton-rich conditions in the wind. 

Neutrino energies depend on the temperature of the medium in the region where neu- 
trinos decouple from matter. This region is known as the neutrinosphere and its location 
is different for each neutrino flavor and energy. Inside of the neutrinosphere, neutrinos are 
in the thermal equilibrium by charged-current reactions (Eqs. (j3J)-flll)). Outside of it, neu- 
trinos escape while their temperature stays almost constant and approximately equal to the 
temperature at their neutrinosphere. Because neutrons are more abundant than protons in 
the neutron star, electron neutrinos continue interacting to larger radii and thus to lower 
temperatures than antineutrinos (gp e > e„ e ). The /x and r (anti) neutrinos interact only via 
neutral-current reactions and decouple at smaller radii, therefore their energies are larger. 
During the first seconds after the explosion, the proto-neutron star deleptonizes and the 
amount of protons in the outer layers decreases. The electron antineutrino energies are thus 
expected to get higher than the electron neutrino energies leading to an increase of the ra- 
tio gp./gj/.. New hyd rodynamical simulations including detailed neutrino transport (e.g. 



Hiidepohl et al.l 120101 ) show that this simple picture is not valid as the spectra of electron 
neutrinos and antineutrinos are rather similar, probably due to neutral-current reactions 
that act in a similar way on all neutrino flavors and b ecome more important as the neutron 
star cools (JHudepohl et al .11201 Ol : lArcones et al.ll2008l ). If electron neutrino and antineutrino 
energies are similar, the neutron-to-proton mass difference favors a value of Y e > 0.5 (as 
shown in Fig. [I]) by reducing the antineutrino absorption rate relative to the neutrino rate. 



-7- 













> 


/ / 








<u 


/ / 




// m 




— 20 








^^^^^ 


>. 


*y 








ai 






0^ 




CD 










c 










<u 15 






+ 


^ -""''^ 


o 










c 




A 


^^ 


_- 


13 










<¥ 10 








-^ol-^ 


4-1 










c 










(T3 











10 12 14 16 

neutrino energy [MeV] 



18 



20 



Fig. 1. — Electron fraction based on Eq. (|6]) for different neutrino (x-axis) and antineu- 
trino (y-axis) energies. The black contours correspond to L De /L Ve = 1 and the grey con- 
tours to L Pe /L Ue = 1.1. The symbols show the electron neutrino and antineutrino en- 
erg ies {e v „ ~ 4.1/ c Xp f or different supernova models re ported in the l iterat ure: square 
for I Woosley et al.l (J19941 ). circle for model M15-ll-r6 of lArcones et al.l ( 120071 ) (note that 



these models do not include red shift corrections which would reduce the neutrino energies 
close to th e neutron staA triable for a 10 M Q progenitor of hsener et J Ml. M H Ji- 



amond for iHiidepohl et al.l (120101 ). all at 10 s after bounce. Adapted from 
Jl996h . 



Qian fc Woosley 



It should be noted that some uncertainty remain in the calculation of the neutrino spec- 
tra. At the high temperatures and densities of the outer layers of the prot o-neutron star, 



there may also be a non-negligible abundance of light nuclei: 2 H, 3 H, and 3 He (lArcones et al. 



2010). These details in the composition are not included in state-of-the-art supernova simu- 
lations, although neutrino in teractions with these light nuclei have an impact on the neutrino 
spectra and thus on the Y e (lArcones et al.l 120081 ). 



2.2. Wind simulations 



In order to accurately calculate the integrated nucleosynthesis in neutrino-driven winds, 
the evolution of the supernova ejecta has to be followed for several seconds using hydrody- 
namical simulations. Such long-term modelin g is currently difficult since the supernova explo- 
sion mechanism is not yet well understood ( IJanka et al.l 120071 : iNordhaus et al.l 120101 ) and it 
is computationally expensive to perform long-time, multidimensional, systematic studies for 
different progenitor stars, as would be desirable in nucleosynthesis studies. Ways to overcome 



these problems include using parametric steady-state wind models (e.g. iThompson et al. 



2003 



2001) and forcing an explo si on by artificially changing neutrinos pro perties (IMesser et al. 



Frohlich et al.ll2006al Jbl; lArcones et al.ll2007t iFischer et al.ll2010l ). The evolution of the 



outflow is rather independent of the details of the explosion mechanism, but depends more 
on the evolution of the neutron star and on the neutrino emission. Therefore, such approx- 
imations are a good basis for nucleosynthesis studies. Although steady-state wind models 
cannot consistently describe hydrodynamical effects (like the r everse shock and mu ltidimen- 
sional instabilities), both approaches agree in the wind phase ( lArcones et al.l 120071 ) . 



For the nucle osynthesis stud i es in this paper, trajectories, i.e. density and temperature 



tonian hydrodynamics (JScheck et al. 



evolutions, from lArcones et al.l (12007) were used. These simul ations are based on New- 



2006 



neutrino transport treatment (JScheck et al. 



corrections for the gravitationa l potential (Marek et al.l 120061 ) combined with a simplified 



Kifonid is et al.l 120061) with general relativistic 



20061 ). which is computationally very efficient 



and reproduces the results of Boltzmann transport simulations qualitatively. The neutrino 
spectra are assumed to follow Fermi-Dirac distributions with spectral temperatures different 
from the local matter temperature in general. The central part (p > 10 13 g/cm ) of the proto- 
neutron star is replaced by a Lagrangian inner boundary placed below the neutrinosphere. 
This reduces the computational time and is justified in part due to the uncertainties in the 
high-density EoS. The evolution of the inner boundary is parametrized by its radius and 
neutrino luminosity. The latter is chosen such that an explosion energy around 10 51 erg is 
obtained. The neutrino luminosity is constant during the first second after bounce and fo l- 
lows a power law decrease afterwards. We use following models from lArcones et al.l (120071 ): 
M10-ll-rl for a 10M o progenitor star (hereafter referred as 10M model), M15-ll-rl (15M 
model) and M15-ll-r6 (15M(s) model) both for a 15M , and M25-15-r4 for a 25M (25M 
model). The proto-neutron star radius and neutrino luminosity evolution are similar for 
models M10, M15, and M25. Model M15(s) has a significant different proto-neutron star 
radius contraction but same neutrino luminosity. While the evolution of the inner boundary 
in model M15 le ads to a more compact neu tron star (approximately following the behaviour 
of the EoS from iLattimer &: Swestyl (Il99ll ) ), model M 1 5(s) r eproduces a neutron star with a 
larger final radius (similar to the EoS from IShen et al.l (119981 )). The different compactness of 
the neutron star implies changes in its radius and mass that directly affect wind parameters 
(Eqs. CP)-©). 

The evolution of radius, entropy, temperature, and density is shown in Fig. [2] for a 
mass element ejected at 5 s after bounce. Since the proto-neutron star in models 10M, 
15M, and 25M follows a very similar contraction and neutrino cooling evolution, S and 
t of these models and therefore their density and temperature evolution are also similar. 
Their wind profiles have a small dependence on the mass progenitor. For the most massive 
progenitor (25M) the proto-neutron star mass becomes larger (due to the higher accretion 



-9- 




10 F 



13 



10 



10" J 



10' 10 

time [s] 







^*^"v. ■> ' 


— 15M 


^O. 


- - 10M 




15M(s) 




-•- 25M 






>^ 1 " ^""^-^i 







10' 



10" 1 
time [s] 



10 u 



400 
350 

^ 300 

c 

3 25 ° 

a. 200 
o 


1 


9, 150 


.. 1 . 


100 


\\ i ■ 


50 




i i i 



10"' 



10" 1 10" 

time [s] 




10' 10" 

time [s] 



Fig. 2. — Evolution of radius, entropy, temperature, and density of a mass element ejected 
5s after bounce for models 10M, 15M, 15M(s), and 25M. 



rate) and consequently the wind entropy is also higher than for the less massive progenitors 
(see Eq. (JT])). The impact of varying the progenitor mass becomes however significant in the 
interaction of the wind with the slowly moving, dense supernova ejecta. This interaction 
results in a wind termination shock or reverse shock. The behavior of the reverse shock 
depends on the wi nd properties, but als o on the progenitor struct ure and anisotropies of the 
early ejecta (see lArcones et al.l (120071 ) : lArcones fc Jankal (120 111 ) for more details). Model 
15M(s), with a slow contracting and less compact proto-neutron star, leads to significantly 
different wind profiles as shown Fig. |2j The expansion timescale is longer and the wind 
entropy lower than model 15M (Eqs. [TH2]). 
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2.3. Nucleosynthesis network 



The evolutions of mass elements ejected between 1 s and 10 s after bounce are used for 
our nucleosynthesis calculations. The final integrated abundances are calculated by adding 
the abundances from each mass element weighted by its ejected mass. Each calculation 
starts when the temperature decreases below 10 GK. The composition is calculated initially 
by assuming nuclear statistical equilibrium (NSE) for a give n initial Y e . The evol ution of 
the composition is followed using a full reaction network ( JFrohlich et al.l l2006bl ). which 
includes 4053 nuclei from H to Hf including both neutron- and proton-rich isotopes. Re- 
actions with neutral a nd charged particles are taken from calculation with the statistical 
code NON-SMOKE R (IRauscher fc Thielemann! 120001 ) and experimental rates are included 



(jAngulo et al 



same as in 



NuDat2| (J2009|). 



1999 , NACRE) whe n available. The theoretical weak interaction rates are the 



Frohlich et al.l ( I2006bl ). Experimental beta-decay rates are used when available 



Results 



The integrated abundances for the wind models presented in Sect. 12.21 are shown in 
Fig. [3j The initial electron fraction and neutrino properties (L u and (e u )) are taken directly 
from the wind simulations. Based on our simulations, elements are produced up to Z = 41 
and no heavy r-process elements ca n be sy nthesized. Figure [3] also shows the LEPP elemental 
pattern from lOian fc Wasserburd (120081 ) normalized to the calculated Z = 39 abundance. 
As the proto-neutron star evolution and wind properties are very similar in models 10M, 
15M, and 25M (see Sect. 12. ip . their final abundances are analogous. The impact of the wind 
parameters (S and r) on the abundances can be seen in the right panel of Fig. [3j Models 
15M and 15M(s) correspond to the same 15 M progenitor, but with different proto-neutron 
star evolution (see Sect. 12. 2p . For model 15M the wind expansion is faster and the entropy 
becomes higher than in model 15M(s). Model 15M(s) can produce mainly iron-group nuclei 
while model 15M leads to heavier nuclei up to Z=47. 

The isotopic production factor indicating the contribution of the ejected isotope to the 
solar system abundances is defined as, 



p® = E 



Mi Xj(t) 



Mf X P 



(7) 



Q\ 



The sum goes over every mass element and Ml- is the total ejected mass obtained after sub- 
tracting the proto-neutron star mass (M ns « 1.4M ) from the progenitor mass. Xj(i) is the 
mass fraction of the isotope % produced in the mass element j, and X Q ^ is the mass fraction 
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Fig. 3. — Integ rated abundances for the models introduced in Sect. 12.21 compared to the 
LEPP pattern (jQian fc Wasserburg||2008l ) rescaled to Z = 39. The abundances of different 
progenitors with a similar evolution of the proto-neutron star are shown in the left panel, 
while the right panel gives the abundances of the same progenitor with different proto- 
neutron star evolutions. 



of the same isotope i in the solar system (JLoddersll2003l ). Following iQian fc Wasserburg 
( 120081 ). the solar system was formed from gas with a contribution of core-collapse supernovae 
of M ras ps 3.3 x 10 10 M Q ejected du ring t ~ 10 10 yr. The r ate of all core-collapse supernovae 



'Kits 



in the Galaxy is .Rsn ~ 10~ 2 yr -1 (jCappellaro et al.lll999l ). This constrains the amount of a 
given element i that can be ejected by a single event. If such an event is representative for 
all supernovae from different progenitors, this implies that 



M 



E^^Ni(M^)<M gas X e 



ej 



©W 



has to be satisfied. Here (M* *) is the average ejected mass, which can be calculated using 
Salpeter initial mass function: 

L m" 2,5 mdm / N 

( M prog> = ^5 7T— ~ !3M . (9) 



j s m 2 - 5 dm 

Therefore, (M*?*) = (M prog ) — M ns ~ 11.5M , assuming that all stars in the range 8M < 
M < 25M contribute. Because the s-process also contributes to the solar system inventory, 
the production factor has to be smaller than the limit given by Eq. (jSJ): 



P(i) < 



M. 



gas 



R S Nt(M^) 



30. 



(10) 



This n umber is an estimate indicative and comparable with the limit of iMathews et al. 



( 119921 ) based on the oxygen production in core-collapse supernovae (P(i) < 10). 
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Production factors are shown in Fig. H] for all models. The dotted horizontal line indi- 
cates the production limit given by Eq. ([TO]) . Such limit assumes that every supernova ejects 
the same amount of matter with the same isotopic composition. If only a subset of super- 
nova (e.g. those from progenitors with masses in a determined range, magnetic fields, etc.) 
eject those isotopes, the limit will move to larger values. Note that only model 15M(s) does 
not lead to significant amounts of heavy nuclei because it is producing mainly iron-group 
nuclei. The overproduction around A = 90 is related with the sudden decrease of the alpha 



separation energy at magic number N = 50 ( s ee e.g.. IWoosley et al.lll994[ IWitti et al.lll994j ; 



Hoffman et al.lll996l ; iFreiburghaus et al.lll999l ; iHoffman et al.ll2008l ). The production factors 



reach higher values for low mass progenitors because the relative contribution of mass ejected 
is higher. 





100 110 120 



Fig. 4. — Production factor given by Eq. ([7]) for the different models shown in Fig. |3j The 
horizontal dotted line marks the overproduction limit. Isotopes for the same element are 
connected by lines. 
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3.1. Impact of the electron fraction 



The electron fraction is extremely sensitive to details of the neutrino interactions and 
transport, and therefore its exact value is expected to be more uncertain than the value of 
the other wind parameters. In this section the impact of Y e on the production of LEPP 
elements is explored. In the following, it is assumed that the temperature of the electron 
antineutrinos is kTp e ^4.5 MeV. For the electron neutrinos, kT Ue is calculated with Eq. ([6]) 
given an initial Y e . Neutrino luminosities are kept constant as a function of time for each mass 
element. For mass element ejected at t = 1 s after bounce, L Ve = Lp e « 25 x 10 51 ergs _1 . 
Subsequent mass elements have constant lumino sities L v {t) = L u {t = ls)/£ 3 / 2 where t is 
the time after core bounce (jAr cones et al.l 120071 ) . Variations of Y e are linked to neutrino 
properties which also affect S and r. In the extreme cases where the change in Y e is rather 
large (i.e., Y e = 0.2 or 0.7), the expansion timescale changes only by a factor of two (which 
does not modify its characteristic timescale of few ms) and the entropy varies by less than 
30% (Eqs. ([I])-©), therefore justifying our approach of keeping S and r as given by the 
simulations while varying Y e . 
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Fig. 5. — Seed nuclei (i.e., nuclei heavier than 4 He), neutron, proton, and alpha abundances 
as a function of Y e based on the mass element ejected at 5s after bounce in model 15M. 
Dotted, solid, and dashed lines correspond to temperatures of 8, 5 and 2 GK, respectively. 



The abundances of neutrons, protons, a— particles, and seed nuclei (A > 4) as a function 
of Y e are shown in Fig. |5] for an expanding mass element ejected 5s after bounce of model 
15M. The abundance at T = 8 GK consists mainly of nucleons and a— particles. At T = 
5 GK, the abundance of seed nuclei has significantly increased but a— particles still dominate. 
Captures of a— particles drive the creation of heavier nuclei until approximately T = 2 — 
3 GK. The composition at this point is very important because the subsequent formation of 
heavier nuclei depends on the free nucleon-to-seed ratio. Three different Y e ranges can be 
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Fig. 6. — Dependence of the abundances of representative elements (Sr, Y, Zr, Cd, Ba and 
Eu) on the electron fraction. These abundances result from a mass element ejected at 5s 
after bounce in model 15M. 



distinguished in Figs. E] and [13 proton-rich (Y e > 0.5), rather neutron-rich (Y e < 0.4), and 
slightly neutron-rich (0.4 < Y e < 0.5) conditions. 

For proton-rich conditions, heavier nuclei are formed by proton and a— capture reactions 
for T > 3 — 5 GK. If Y e is large, heavier elements can be synthesized for lower temperatures 
due to the increase of proton and neutron abundances. The amount of neutrons increases 
with Y e , i.e. with the number of free protons (see blue dashed line in Fig. [5]). When Y e > 0.5, 
there are always more protons available than what can be captured, and only very few 
neutrons produced by antineutrino absorption on protons. These neutrons are immediately 
captured allowing more matter to bypass long-liv ed isotopes (bottlenec ks ) such as 64 Ge by 



n, p) reaction s in what is known as the z/p-process (IFrohlich et al.ll2006al lbl; iPruet et al.ll2006 



Wanajol 120061 ). As the temperature decreases below « 3 GK a combination of (n,p), (n, 7), 
and (p, 7) reactions carry the flow to heavier nuclei. After charged-particle reactions end, 
the produced matter decays back to stability, mainly to neutron-deficient isotopes. Only for 
early mass elements, the reverse schock becomes important since its temperature is still in a 
relevant range (T « 2GK) for the z^p-process to occur. The temperature jump at the reverse 
shock favors thus the production of heavier nuclei. The impact of the reverse shock o n the 
z/p-process will be addressed in future work (IWanajo et al.ll2010at iRoberts et al.ll2010l ). The 
heaviest nucleus that can be reached, depends mainly on the specific ne utrino luminosities 



and energies, and in less scale, on entropy and expansion timescale (see IPruet et al.l 12006 



Wanajol 120061 ; iMartinez-Pinedo et al.l 120061 ). but as shown in Fig. El LEPP elements can be 
produced under proton-rich conditions. As we discussed later, abundances in proton-rich 
winds are dominated by p-nuclei. 
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Extreme neutron-rich conditions are naturally more favorable to reach heavy r-nuclei. 
Although these conditions have not been reproduced in standard supernova ejecta, it may 
still be possible to obtain them in scenarios such a s exp losions with high rotati on and mag 



netic fields where a jet for ms (jCameronl l200ll . 120031 ). or in quark nova ( jOuyed et al. 



20021 ; iJaikumar et al.l 120071 ) where a phase transition in the neutron star leads to a di- 



rect ejection of very neut ron-rich matter. In neutron-rich conditions, the reac tion sequence 
4 He(an,7) 9 Be(a,n) 12 C (IWooslev fc Hoffman! 1 19921 ; iFreiburghaus et al.l Il999l ) followed by 
a— captures quickly increases the amount of seed nuclei to amounts larger than in proton-rich 
conditions as shown in Fig. [5j Once charged-particle freeze-out occurs, nuclei are still driven 
by neutron-captures and /3-decays before the density gets too low and the created matter 
/3-decays towards stability. 

When the ejecta are slightly neutron-rich, the 4 ife(an,7) sequence starts to compete 
with the triple alpha-reaction and it dominates for Ye < 0.45. Subsequent a-captures reac- 
tions drive the flow following a path close to or at t he valley of beta-stability ( IWooslev fc Hoffman 
19921 : Iwitti et al.lll994J IFreiburghaus et al.lll999l). At closed N = 50 neutron shell, the al- 
pha separation energy decreases (jMoller et al.lll997l ) preventing a-captures to continue with 
the f ormation of heavier nuclei . This leads to an accumulation of matter around 88 Sr, 89 Y, 
90 Zr JWooslev fc Hoffman! [l992l ; Iwitti et al.lll994l : IFreiburghaus et al.lll999l ). When the tem- 



perature drops below « 3 GK, the Coulomb barrier produces the end of a-reactions and 
the material captures neutrons if the neutron-to-seed ratio is still high enough. When the 
matter is symmetric (Y e = 0.5) only elements up to the iron-group are formed. 



3.2. 



LEPP and Y e evolution 



For the evol ution of the initial electron fraction, we use a parametric approach (similar 



Wanajol 120071 ) with the aim of fitting the abundances observed in LEPP-enriched UMP 



to 

stars. In the following, model 15M is used unless otherwise specified. 



3.2.1. Proton-rich winds 



The integrated abundances corresponding to three different Y e evolutions ar e show n in 
Fig. These evolutions are motivated by the recent results of iHiidepohl et al.l (120 101 ) for 
the explosion and wind evolution of a 8.8 M & stellar progenitor that shows an increasing Y e 
as a function of time. The abundance pattern in the LEPP region 38 < Z < 47 changes only 
slightly for the different cases and fits observations within uncertainties. The robustness of 
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the abundance pattern can be attributed to the consistency in the (p, 7) and (n,p) reaction 
path for the mass elements responsable for the abundances in the LEPP region. Most of those 
elements are produced when Y e > 0.6. Such robustness provides a natural expla nation for the 
consis tency of the suggested UMP star abundance LEPP pattern mentioned in iMontes et al. 
( 120071 ). Figure |8] shows the similarity of the abundance pattern for Y e = 0.65 for two different 
mass elements ejected at early and late times. The only significant difference is the amount 
of mass ejected in each case. 




Fig. 7. — Integrated abundances in model 15M using the three different Y e evolutions shown 
in the inset. Same colors and line styles are used for the abundances and Y e evolutions. The 
LEPP pattern has been rescaled to fit the solid line abundances. 

The main contribution to the integrated abundances come from mass elements ejected 
during the first seconds after bounce due to their large mass outflow. Late mass elements 
eject less mass but produce heavier elements because of the higher Y e . These trends are 
shown in Fig. [8] by the mass-weighted abundances at 1 s and 9 s after bounce for two 
different initial Y e . For Y e = 0.5 mainly iron-group nuclei are produced. LEPP elements 
can be synthesized for Y e = 0.65 as antineutrino absorption on free protons make possible 
(n, p) reactions. As high Y e values are reached only by mass elements emitted at late times 
(small inset in Fig. [7]), LEPP abundances are weighted by a lower mass ejection of the shell. 
Depending on the Y e of the early mass elements, the amount of iron-group elements can 
change by orders of magnitude leading to large variations of the ratio between iron-group 
and LEPP elements. 

While the LEPP elements are most probably produced during the neutrino-driven wind, 
the synthesis of iron-group nuclei occurs mainly in explosive nucl eosynthesis when the super- 
nova shock disrupts the stellar envelope (jWoosley fc Hegerll2007l ). Therefore, the production 
of the iron-group elements in the neutrino-driven wind should be negligible. Although in our 
simulations iron-group elements are produced mostly when Y e « 0.5, a higher Y e combined 



-17- 



io" 3 

10" 6 

10- 7 

10" 8 
10' 9 



E 10" : 
< 10" J 

io- ! 







10 15 20 25 



30 35 

Z 



40 45 



50 



Fig. 8. — Abundances for the 15M trajectories ejected at 1 s (solid lines) and 9 s (dashed 
lines) after bounce assuming that the initial electron fraction is Y e = 0.5 (black lines) and 
Y e = 0.65 (grey lines). 



with a lower neutrino luminosity may also produce them. In our model 15M, the abun- 
dances of ir on-group elements (in particular, Z=29) compared to observed abundances of 
HD 122563 ( jWestin et al.ll2000l ) can constrain the Y e of the early trajectories. We find that 
the Y e evolution has to follow either a rapid increase so the amount of iron-group elements 
ejected stays low, or it has to start already at Y e > 0.55 so there is only a minor production 
of these elements. The evolution represented by the solid line in Fig. [7] produces too much 
iron-group nuclei relative to LEPP elements as the Y e stays close to 0.5 for several seconds. 
The evolutions represented by dotted and dashed lines are consistent with observations as 
the abundances of the iron-group elements are less than 10% of the observed abundances 
when normalizing to Z=39. 

Production factors corresponding to different proton-rich evolutions are shown in Fig. |9j 
To gain some insight into the absolute abundances of LEPP elements, one may compare their 
abundances with the measured values of a-elements that are the dominant products of core- 
collapse supernovae. For instance, Mg should take the value given by Eq. ( fit)]) since it is a 
typical super nova product. Assuming HD 122563 is representative for all supernovae and 
[Mg/Y]=0.9 (IHonda et al.l 12004 120061 ). one would expect the production factor of Y (single 
isotope) to be around 3. Figure M shows that P(Y) has roughly the right value. Nevertheless, 
that some isotopes are over the production limit given by Eq. (Tl0|) means that either not 
all supernovae in the range 8M < M < 25M Q reach the same proton-rich conditions or 
that the ejecta from our spherically symmetric, parametric model cannot be extrapolated to 
every supernova ejecta. If the abundances observed in HD 122563 are not representative of 
a typical supernovae ejecta then Y (or other LEPP elements) could be produced in different 
amounts. We find that abundances from models 10M and 25M behave very similar to 15M. 
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Fig. 9. — Production factors corresponding to the elemental abundances presented in Fig. [7] 
by the dashed and dotted lines are shown in the left and right panels, respectively. Almost 
only neutron-deficient isotopes (p-nuclei) are synthesized. The i sotopes marked with circles 
are not produced in enough quantities in iTravaglio et all (J2004f ). 



For model 15 M(s) is not possib le to produce LEPP elements even when Y e ~ 0.65 due to the 
low entropy (jPruet et al.ll2006l ). In addition, most of the abundances are on the proton-rich 
side of the valley of beta-stability (neutron-defici ent isotopes). If t he LE PP were responsible 



for the isotopic underabundances reported by ITravaglio et all (120041 ). proton- rich winds 



would not be enough to account for this contribution. 

The nucleosynthesis dependance on neutrino luminosity is shown in Fig. [10] for a tra- 
jectory ejected 1 s after bounce with Y e = 0.65. The neutrino luminosity obtained in 



Arcones et al.l (120071 ) has been reduced by different factors. When the neutrino luminos- 
ity decreases the wind entropy has to increase (Eq. [TJ. For our test, we do not consider this 
dependence and keep the entropy (and other wind parameters) unchanged. Although, this 
is not fully consistent, it allows us to investigate the influence of the neutrino luminosity 
on the abundances separately from the effect of other wind parameters. We find that a re- 
duction of the luminosity by a factor 2 (dashed line, "L/2" ) leads to the production of only 
Z < 44 elements. A factor of 5 reduction is enough to inhibit the production of any LEPP 
elements, as the flow can barely pass through the bottlene ck at 64 Ge. We note t hat previous 
nucleosynthesis studies based on the explosion model of Kitaura et al.l (120061 ) do not pro 



2009( 1 because the (anti)neutrino 



duce LEPP elements (jHoffman et al.ll2008l ; IWanajo et al. 
luminosity of this low mass progenitor (M = 8.8M ) is rather small and the expansion is 
very fast. Therefore, there are not enough neutrinos during sufficient time to get a successful 
^p-proc ess. The neutrino luminosities obtained for heavier mass progenitors, as in our case 
and in iPruet et al.l ( 120061 ) for a 15 M & progenitor, show that with higher luminosities the 
z/p-process can occur in the supernova outflows. Although we have shown that with the 
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calculated (anti)neutrino luminosities it is possible to obtain the LEPP pattern, the neu- 
trino luminosities from hydrodynamical simulations have large uncertainty and prevents us 
to draw definite conclusions. 
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Fig. 10. — Abundance patterns obtained for the mass element ejected 1 s after bounce with 
Y e = 0.65 and different neutr i no lu minosity. Solid line shows the abundances using the 
luminosity from lArcones et all ( 120071 ) . For the other cases the luminosity has been reduced 
by different factors as indicated in the label. 



3.2.2. Neutron-rich winds 



Figure [IT] shows integrated abundances for three possible neutron-rich Y e evolutions 
(inset of the same Figure). When the Y e decreases from 0.5 to 0.486, only Z < 40 are 
produced (similar to Fig. [3]). If Y e decreases to 0.45, elements up to Z = 47 can be produced. 
Contrary to proton-rich conditions, different Y e evolutions lead to significant changes in the 
pattern. This sensitivity in the final abundances is due to the abrupt drop of the neutron 
abundance near Y e = 0.5 (shown in Fig. [5]). This rapid change in neutron abundance is 
accompanied by a non-linear change in LEPP abundances when 0.4 < Y e < 0.5 (shown in 
Fig. [6]). Moreover, mass elements ejected at slightly different times (and therefore with small 
variations of the entropy and the expansion timescales) show a substantial difference on the 
neutron-to-seed ratio and subsequently on the abundances, making again difficult to get a 
robust pattern. Such strong dependence of the final abundances on the wind parameters is 
also found for the other models introduced in Sect. El In model 15M(s) the electron fraction 
has to be reduced down to Y e w 0.4 to produce elements up to Z = 47. 

Production factors of the abundances shown by dashed and dotted lines in Fig. [11] are 
presented in the left and right panels of Fig. [121 respectively. Since in neutron-rich winds 
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Fig. 11. — Same as Fig. [7] but for neutron-rich Y e evolutions. 
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Fig. 12. — Production factors for two evolutions of the electron fraction: left panel correspond 
to green dashed line and right panel to the red dotted line in Fig. [HI The dotted horizontal 
line represents the upper limit above which isotopes are overproduced (see Sect. 13.11) . 



(0.45 < Y e < 0.5) the amount of heavy nuclei produced is higher than in the proton-rich 
winds (Fi g. El), the overproduc tion around A = 90 (related to the neutron magic number 



N = 50, iHoffman et al. 



(119961 )) becomes a constraint. Similar to the conclusion reached 



for the proton-rich conditions, this suggests that not all mass ejected in neutrino-driven 
win ds in core-collapse supernoyae experience neutron-rich conditions. Previous results (s ee 



.g.. IWooslev fc Hoffman 



Freiburghaus et al.lll999 : 



1992 



Woosley et a 



19941 ; IWitti et alJll994J ; IHoffman et al.lll996 



Hoffman et al.ll2008l ) suggest that only a small amount of the mass 



ejected by core-collapse supernovae can be neutron-rich. It is possible that the neutron 
richness of the ejecta depends on progenitor mass because more massive stars lead to more 
compact neutron stars and thus higher neutrino energies, which may favour neutron-rich con- 
ditions. A different alternative to avoid the overproduction is to have anisotropic explosions 
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with an electron fraction distribution in the ejecta that achieves neutron-rich conditions only 
in some regions. Recent two-dimensional supernova simulations pr esent proton- rich ejec ta 
with neutron-rich blobs or pockets that contain only a small mass ( IWanajo et al.l 1201 Obi ). 



Further investigations are necessary to conclude which of the possibilities addressed here 
can explain the observations of LEPP elements in UMP stars. Based on spherical symmetric 
simulations we have showed that it is possible to reproduce the LEPP elemental pattern 
under different wind conditions. Multidimensional simulations will determine the effect of 
anisotropic Y e distribution on the wind nucleosynthesis. Such multidimensional simulations 
combined with exploratory studies like ours will provide constraints for the specific conditions 
(e.g., Y e evolution) and for the overall contribution of neutron- and proton-rich winds to 
the solar system inventory. Using our re sults, it is hard to obtain the solar LEPP pattern 
(JTravaglio et al.ll2004j ; iMontes et al.ll20071 ) by combining neutron- and proton-rich conditions. 
Moreover, Galactic chemical evolution models are necessary before a definite conclusion 
can be drawn to the exact contribution of core-collapse supernovae to the solar systems 
abundances. 



4. Conclusions and summary 

We have perf ormed nucleosynthes is calculations based on spherically symmetric super- 
nova simulations (jArcones et al.ll2007l ) which include efficient neutrino transport and follow 
the evolution of the ejecta from a few milliseconds to several seconds after bounce during 
the neutrino-driven wind phase. The integrated nucleosynthesis was analyzed for different 
stellar progenitors with masses of 10, 15, and 25 M & . We have found that light element 
primary process abundances can be produced under realistic conditions in the neutrino- wind 
phase. No heavy r-process elements can be synthesized under those conditions (see also 
Roberts et allboiol : IWanaio et aJboiObh . 



Our results indicate that the nucleosynthesis occurring in neutrino-driven winds does not 
depend significantly on the progenitor, but rather on the proto-neutr on star evolution which 
in turn directly affects nucleosynthesis-relevant wind parameters (jQian fc Woosleyl Il996l ) 
such as entropy, expansion timescale, and electron fraction. The electron fraction depends 
on neutrino properties which are determined by the still uncertain neutrino interactions and 



trans port. Since the calculation of the electron fraction remains challenging (jHudepohl et al. 



20101 ) the impact of the electron f raction on the produ c tion o f LEPP elemen t s was studied. 
The LEPP pattern obtained by iQian &: Wasserburgl (120071 ); IMontes et al.l (120071 ) was re- 
produced for different evolutions of the electron fraction towards proton- and neutron-rich 
conditions. 
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The LEPP pattern reproduced in proton-rich winds was found to be robust under small 
variations in the evolution wind parameters. Nevertheless, elements heavier than iron-group 
nuclei can only be produced when the neutrino fluxes are high enough to allow for a success- 
ful ^p-process. The amount of heavy elements ejected is rather low since most of the outflow 
matter consists of protons and a-particles. Although the elemental abundances nicely repro- 
duce the observed LE P P pattern, mainly ne utron-deficient isotopes are produced (see also 
Frohlich et al.ll2006al Jbl; IWanajo et al.ll2010al ). Therefore, proton-rich conditions can explain 
the LEPP elements observed in UMP star s but not in the missing isotopic underabundances 
in the solar system (JTravaglio et al.ll2004j ). 



When the electron fraction is assumed to evolve towards neutron-rich conditions, the 
LEPP pattern can be also reproduced but it is not robust under small variations of the wind 
parameters. We find an overpro duction around A ~ 90 that was already po i nted out in 
previous nucleosynthesis studies (IWooslev fc Hoffmanlll992l ; IWoosley et al.lll994j ; IWitti et al. 



1994] ; iFreiburghaus et al.lll999t iHoffman et al. 



2008). T his overproduction an d the fact that 



most recent supernova simulations (IFischer et al.ll2010l ; iHiidepohl et al.ll2010l ) favor proton- 
rich winds could suggest that neutron-rich winds are rare events. Scenarios that could explain 
the ejection of only a small amount of neutron-rich material include 1) very small pockets 
with n eutron-rich material that may appear in multidimensional simulations (IWanajo et al. 
2010b(); 2) neutron-rich conditions that can only be obtained at very late times when the 
amount of ejected mass is very small; 3) there is only a small subset of all supernovae that 
can develop neutron-rich winds. 



The isotopic underabundances reported by (JTravaglio et al.ll2004l ) cannot be reproduced 
based in our models without overproducing other isotopes. Evolving from proton to neutron 
conditions in the wind does not resolve the overproduction, which could be a hint of the 
limitatio n of our models. How ever, this can also suggest deficiencies in the s-process model 
used by (JTravaglio et al.ll2004l ). We cannot conclude whether the solar and the UMP stars 
LEPP abundances are produced by the same light element primary process. 

Observation of isotopic abundances in UMP stars would constrain the evolution of the 
electron fraction and thus of the neutrino properties in supernovae. In addition more work 
is necessary in improving long-time supernova models (e.g., multidimensional simulations, 
more accurate treatment of neutrino reactions and transport) and also the experimental and 
theoretical nuclear reactions relevant in the wind nucleosynthesis before definite conclusions 
can be reached. 
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